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We prove the multi-band Bloch oscillation and Stark ladder in the nk and site representation from 
the Floquet theorem. The proof is also possible from the equivalence between the Floquet system, 
Bloch oscillation, and the rotator with spin. We also exactly solve the periodically driven two level 
atom and two band Bloch oscillation in terms of Heun function. 
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I. INTRODUCTION 

Periodically driven two level atom with Rabi oscillation 
effect is one of the simplest quantum systems, yet it still 
has not been exactly solved [l|, 0, d, 0, H, H, Q • The semi- 
classical Hamiltonian is the interaction between a time 
periodic field and a two level atom. It is the approx- 
imation to the interaction between a quantized photon 
field and a two level atom when the the photon num- 
ber is large Q. Floquet system (FS) has a time periodic 
Hamiltonian. Periodically driven two level atom is the 
simplest FS. 

Bloch oscillation (BO) and Stark ladder (SL ) are an 
old and controversial problem [1, i, [l^, [11 Q 

d, [M [13, H [H, H [In . One band BO and SL are easy 
to be established, but as shown by Zak [13], there is a 
paradox in the one band BO. Other proofs of BO and 
SL in the multi-band cases or in the case of the general 
BO Hamiltonian includes [3, but one proof 

is general followed by some comment and reply. Nen- 
ciu argued when taking into consideration the inter-band 
hopping, SL does not exist any more [1^. In the paper, 
we establish the multi-band BO and SL from the Flo- 
quet theorem. Even if our proof is not general enough to 
end the controversy, at least, it points out the physical 
mechanism of BO and SL: the Floquet theorem. 

The BO Hamiltonian in the nk representation allows 
an interpretation as a FS. Based on the Floquet theorem, 
we give a proof of BO and SL both in the nk and site 
representation. 



II. EXACT SOLUTION OF A PERIODICALLY 
DRIVEN TWO LEVEL ATOM 

The Hamiltonian of a periodically driven two level 
atom is 



H{t) = ri(T^ — Aux sinwt, 



(1) 



where Gx and are Pauli matrices, and and A are 
parameters. After a time-independent unitary rotation 
in the Hilbert space 22], the Hamiltonian becomes 



The Schrodinger equation is 



[Auz sinwi -I- VIgx) 



(3) 



where (ci, 02)"^ is the wave function of the two level atom 
after the unitary transform. For simplicity, we set h = 1 
and uj = 1. 



d 
d 

[i— + As\TLt)c2 = rici. 
at 



After removing C2, we get 
92 



9t2 



ci -f {lA cos t + A'' sin^ t + Q'')ci = 0. 



(4) 



(5) 



We change the variable from t to z. z = ^ cost + i. We 
get 



(J2 ^ g 

~ '^d^'' + - 2^0;"' 

+ {AA^z^ - AA^z - 2iAz + iA- n^)ci = 0. 

The general solution is 

ci{z)=e^''^'x 

Hc{AiA,-\,-\,~2iA,iA+l - VL^,z) 



(6) 



(7) 



^/zHci-iiA, \,-\,-2iA, iA 



n^z) 



where He is the Heun confluent function. See the ap- 
pendix of the Heun function. The general solution is 

ci(0 = e2"-4(cost+i)x 

Hc{AiA,-\,-^,-2iA,iA+l - n'^,\cost+ \) ('g^ 
jyc(4iA, i,-i,-2iA,iA-h I - f]2 1 cost + i)x 



(cost 1)1/2, 

When t = 0, z = 1, which is a singular point of He. Now 
we do not have enough knowledge to handle the singular 
point in order to consider the initial conditions. When 



H'{t) = AcTz sin Ltit + flcTx- 



(2) 



H'{t) = A(Tx cost + ilaz 



(9) 
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The general solution is 



The eigenstates are 



"^^^^^ ^ Hd(-4A, -41^2, -8A,4f]2,i cot |)e^^"'"* 



(10) 



III. EXACT SOLUTION OF TWO BAND 
BLOCH OSCILLATION 



Ank~it sin A; 



i(A± VA2 + (ef^i?)2) 



(17) 



The result is first derived by Fukuyama et al in 14]. The 
physical meaning of Eq. (16, 17) is two SLs. Electron 
oscillates between two bands with a period -^^==5. 



The Hamiltonian of two band Bloch electron in a elec- 
tric field is 



B. Exact solution of two band Bloch oscillation 



H= J2 [{^-eFn)aian + {-A-eFn)bibn 

n— — 00 



3FR{albn + b^an)], 



(11) 



When [Hi^, Hs] 7^ 0, the rotation and spin degrees of 
freedom are coupled together. The Hamiltonian eigen- 
value equation is 



HiP = Ei;. 



(18) 



where a and b refer to electrons in two bands with band- 
widths 2ta and 2tf, respectively; the first two terms are 
the site energies, the middle two describe site-to-site hop- 
ping, and the last term gives interband hopping Il4l. 12111 . 

The Hamiltonian can be Fourier transformed [ij, l23| 
into the nk representation [l^ lUi HE]) in which the 
Hamiltonian is 



H = 



A - ieF 



dk 

eFR 



ta cos k —eFR 

—A — ieF-g^ + tft cos k 



[12) 



The eigenstates of the Hamiltonian is periodic of k \\i 
We assume the period is 2-k. 
We define 



ta COS k 

tb COS k 



as the rotation Hamiltonian and 



A -eFR 
eFR -A 



(13) 



(14) 



as the spin Hamiltonian. The names will be explained in 
the Sec. IV. 



A - ieF 



-§Z+ta cos k 

eFR -A-ieF-^ + tbCOsk 



-eFR 

dk 



E 



(19) 



where ijj = (ci,C2)"^ is the wave function and ci and C2 
are functions of k. 

d 

(-ieF— +taCosk + A- E)ci - eFRc2 = 0, 
ok 

d 

-eFRci + {-ieF—+tbC0sk-A-E)c2 = 0.(20) 
We set the unit eF = 1. After removing C2, 



92 d 
g^ci + {{ta + tb) cos k - '2E)i—ci 

ita sin k — tatb cos^ k — (ia(— A — E) 
- tb{A - E)) cos k-E"^ + A^ + rA ci = 0. 



(21) 



A. Rotation and spin decoupled 

Now we discuss the simplest two band BO. If ta = ^6 



d 

H = —ieF— + t cos k - 
ok 



A -eFR 
-eFR -A 



(15) 



Then [Hr, Hs] — 0. The rotation and spin degrees of 
freedom are decoupled. The eigenvalues of Eq. (15) is 



^± = neF ± + (eFi?)2. 



(16) 



We change the variable from ci(fc) to u{k) and define 
ta^tb — S. ci{k) — u{k) exp(—iasmk). 



Q2 q 

TTT^u - (S COS k + 2E)i—u 



dk 

5{A - E) cosk + A'^ + R^ - E'^ ]u ^ 



(22) 



Eq. (22) can be solve by Heun function. What we have 
to do is to transform it into the standard form of Heun 
function. We change variable from fc to z. z — / cot ^. 
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+ [4(i;2 - - i?2 _ + ^E)z'^ + Ai-E"^ + + i?2 - + 5E)]u = 0. 
After another variable change, Eq. (23) can be written into the standard form. The general solution is 



l + z e \ Hd (2(5, -4(5 A - 4A2 - 4i?2 , -4(5, -4(5A + 4A2 + 4i?' 



z) 



""^^^ " ^ 1 + z ^ \Hd (-2(5, -4(5A - 4A2 - m\ -4(5, -4(5A + 4A2 + 4i?2, z) e-25./(.^-i) J- (24) 

^42<^'-4'5A-4A2-4i?2,-45,-45A + 4A2 + 4i?2,icot(f)) \ 

-e i74-2(5,-4(5A-4A2-4i?2,-4(5,-4(5A + 4A2 + 4i?^icot(|))e''5™'= /■ ^ ' 

(h\- »bJ (2^, -45A-4A2-4i?2, -4(5, -4(5A + 4A2 + 4i?2, zcot(|))e-^*» 1 , . 

^^•^ \ i/4-2(5,-4(5A-4A2-4i?2,-4(5,-4(5A + 4A2 + 4i?2,zcot(f))e-'*''^"'= /• ^ ^ 



IV. FLOQUET SYSTEM, BLOCH 
OSCILLATION, AND STARK LADDER 

The Schrodinger equation of a FS is 



(27) 



where 5'(t) is the wave function and H{t) = H{t + 2tt). 
Due to the Floquet theorem, the fundamental solutions 
of Eq. (27) are the multiple of a c-number and a time 
periodic wave function. 



(28) 



where is a real number and referred as quasienergy Q . 
It can be verified 



(Hit) - ^h^^-^Mt) = Em- 



(29) 



Now we change the variable from the time t to the Bloch 
number k. t = k. 



(H(k) - ih-^mk) = E<^{k). 

ok 



(30) 



$(fc) is seen as the wave function of the Bloch electron. 
The general BO Hamilton 



H ^ — + Vix) - eFx, 
2m 



(31) 



where m,x,p are the mass, position and momentum of 
the Bloch electron. In the nk representatioUja multi- 
band approximation of Eq. (31) is [H, [H, [M Isl 



M 



M 



d 



ff(fc) = ^e„(fc)- eFXn,n{k)-ieF—, (32) 



n.m— 1 



where X„„i(fc) 



(33) 



is periodic of k. 

Eq. (30) and (32) have the same form. The eigenvalue 
problem of the BO Hamiltonian is 



M 



M 



^e„(fc)- ^ eFXnrnik) - ieF 



n=l 



n,m— 1 



d_ 

dk 



$(fc) = E^{k), 
(34) 

If {E,^{k)} is the eigenvalue and eigenstate of Eq. 
(30) and (32), then 

E' = E + m X 2TTEe, (35) 

$'(fc) = e*'"^2^'^$(fc), (36) 

gives another eigenvalue and eigenstate of the BO Hamil- 
tonian. But it gives the same solution to the FS Eq. (27). 
In FS, if H{t) is a M X M matrix, which corresponds 
to a M band BO Hamiltonian, the total number of the 
quasienergy is M. The eigenvalues of the BO Hamilto- 
nian are grouped into M SLs. If M is finite, BO can not 
have continuous spectra and the eigenstates of BO Hamil- 
tonian are localized. If the electron is put on one site, 
the wave function oscillates because the wave function 
can be expanded by (approximately) finite eigenstates of 
the BO Hamiltonian. It is the connection of Eq. (32) 
to the Floquet system, that gives the eigenvalues of SLs. 
We note Fukuyama et al first used the Floquet theorem 
to prove SL in a two band case Avron et al also 

gave a proof of SL based on Eq. (32) in 15]. But our 
proof is clearer. 

We require H{k) is not a function of i^. In this way, 
we implicitly assume 



u*j^k{x)umkix) dx = 



(37) 
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or 

(ns|ms) = (38) 

when n ^ m in the site representation. Under this 
assumption, the interband hopping matrix elements 
Xnmik) or 

{ns\ — eF{x)\ms) — {ns\ — eF{x + Xo)\ms) (39) 
do not depends on x. 

A. Infinite matrix representation 

Every FS and BO Hamiltonian can be represented as 
an infinite matrix. The infinite matrix of the periodically 



driven two level atom was first written out and referred 
as Floquet Hamiltonian by Shirley The FS formal- 
ism [3,13 including the infinite matrix representation is 
widely used in the practical numerical calculations, such 
as atomic and molecular excitation, ionization in a laser 
field The Hamiltonian of a periodically driven two 
level atom is 

H{t) = Vla^ + Acos{Ljjt)a^. (40) 
Then the Floquet Hamiltonian 

HF = H{t)-i^^, (41) 

has the following infinite matrix representation 0, 
in the \sn) = \ t, 7=6-'"*) or | i, ^e"™*) basis 





A/2 




















A/2 


Q + u 








A/2 



















A/2 




















A/2 


n 








A/2 








A/2 










A/2 




















A/2 


n — (jj 

















A/2 








-n - 


-Lj A/2 




















A/2 


Q.-2UJ 



(42) 



which has a more compact form 



Hs + 2uj 

Hr 









Hr 

Hs+uj 

Hr 








Hr 

Hs 

Hr 








Hr 

Hs' 

Hr 







Hr 

Hs- 



2uj 



where 



n 

-n 



A/2 
A/2 



(43) 
(44) 



The infinite matrix representation of the two band BO 
Hamiltonian Eq. (11) is 



Hs + 2eF 
Hr 






Hr 

Hs + eF 

Hr 








Hr 

Hs 

Hr 








Hr 

Hs- 

Hr 



eF 







Hr 

Hs ' 



2eF 



where 



A 

-eFR 



-eFR 
-A 



; Hb — 



ta/2 







4/2 



(45) 



(46) 



Note in this section Hf> is different from the Sec. HI by 
the part of —ieF-^. 

The Floquet Hamiltonian of periodically driven two 
level atom Eq. (43, 44) and the two band BO Hamilto- 
nian Eq. (45,46) are very similar. The difference is just 
the different basis of the Hilbert space and parameters. 
Given a FS, we can find a BO Hamiltonian, which has 
the same infinite matrix representation as the Floquet 
Hamiltonian of the FS and vice versa. FS must has very 
similar dynamic behaviors with BO and vice versa. This 
is the reason why we can use the Floquet theorem to 
prove BO. 



B. Floquet system, Bloch oscillation, and rotator 
with spin 

The Hamiltonian Eq. (32) is an approximation of Eq. 
(31) based on the Bloch function. The proof of BO and 
SL is achieved in the nk representation. But Eq. (31) 
can also be approximated in the site representation just 
as two band approximation Eq. (11). 

The BO Hamiltonian of a multi-band Bloch electron 
in a linear electric field is a rotator with the spin degree 
of freedom (RS). The basis of the Hilbert space of Eq. 
(11) is \sn), where s is one of the multi-bands and n is 
n-th site. In the rotator representation, the basis of the 
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Hilbert space is \sn), where s is the spin of the rotator and 
n is i^e"'"*^ with k the angle of the rotator. In the Sec. 
III, we rewrite Eq. (11) as Eq. (12) from the perspective 
of the Fourier transform. But this can also be achieved 
by the site-band and rotation-spin correspondence. The 
n-th site corresponds to -^=e~'^"'^ and the s-th band the 
S-th spin state. From this correspondence, we can rewrite 
Eq. (11) into Eq. (12) without the Fourier transform. 
Every BO Hamiltonian can be rewritten as a RS. H{k), 
such as in the Eq. (12), is the coupling between rotation- 
rotation, spin-spin and rotation-spin degrees of freedom. 
The Bloch wave number in the nk representation k is 
the angle of the rotator. By corresponding the Hilbert 
space basis and the Hamiltonian, we have established the 
equivalence between BO and RS. 

For the Hamiltonian of every one of FS, BO and RS, 
we can find another two Hamiltonians in the other two 
systems, and the former Hamiltonian has the equivalent 
behavior with the latter two. 

Generalizing Eq. (12), in the rotator representation, 
Eq. (31) is approximated as 

M 



M oo Q 

- ^ ^ ^ eFXsn,s'm\sn){s''m\~ ieF—, 

s,s' — l n,m— — GO sn^s'm 

(47) 

where \sn) — \s, —^e~'^'^^) is the basis of the RS and 



Xsn,s'm^ I U*nix)xUs'm{x)dx (48) 



with Uns and u^s being Wannier functions, give the cou- 
pling between bands and sites. The specific value of 
Xsn,s'm is gave by the parameters of the material and 
the electric field, but it is "simulated" by the RS. We 
note the rotator with a linear — i-^ in the Hamiltonian 
is first studied by Grempel et al [23| and Berry [28|. 

We assume Xsn.s'm — -^sn+n' ^s'm-\-n' ^sO.s'm— n- 



of periodically driven two level atom and two band BO 
in terms of Heun function. 

The fully quantized field and two level atom interac- 
tion is more interesting. But the present method com- 
pletely fails because first, after the second quantization, 
the Hamiltonian is not doubly infinite; second, the off- 
diagonal matrix changes away from the ground state . 

The original BO Hamiltonian Eq. (31) can be ap- 
proximated in the nk representation as in Eq. (32) or 
in the site representation as in Eq. (11). In the nk 
representation, we can prove BO and SL from the Flo- 
quet theorem. In the site representation, we first rewrite 
the same Hamiltonian in the rotator representation, then 
again prove BO from the Floquet theorem. We do not 
know whether the the two methods are equivalent. Al- 
though, the nk representation can be Fourier transformed 
into the site representation, it is difficult to establish the 
equivalence between Eq. (33) and (48). Both approxi- 
mations may have their own merits in applications. We 
note the multi-band approximation seems impossible to 
remove the paradox of Zak [13] . 

Our method can not treat the most general BO prob- 
lem with the Hamiltonian Eq. (31). But in a physi- 
cal material, the Bloch electron only occupies the lowest 
several bands. So we believe the BO and SL are highly 
possible in reality. 

Our most significant contribution in the paper is the 
equivalence between the FS, BO and RS. We think the 
three are the same problem with three "faces". The 
equivalence between FS and BO is established by com- 
paring the Floquet Hamiltonian Eq. (29) and the BO 
Hamiltonian Eq. (30,32). The equivalence between BO 
and RS is established by corresponding the site-band ba- 
sis with the rotation-spin basis. At last, the three share 
the same infinite matrix structure. 

The equivalence is shown as the following figure. 



M 



i7^^e,(fc)|s)(a 



s=l 
M 



n— — oo 



dk' 



(49) 

From Eq. (49), we can also prove BO and SL based on 
the Floquet theorem. 




CONCLUSION AND DISCUSSION 



In summary, we have given a proof of BO and SL based 
on the Floquet theorem. We also give the exact solution 
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APPENDIX A: HEUN FUNCTION 

The Heun function H{a, q; a, (3, 7, S, z) is defined as the 
solution of the following equations 



7 



53 + <T 



— q 



a dz z{z — l)(z — a) 



H{a,q;a,f3,j,S, z) = 0, 
H{a,q;a,f3,j,S,0) = 1, 
dH{a, q; a, /3, 7, (5, z) 



dz 



|z=0 



07 ' 



(Al) 



where e = a + /3 + l — 7 — 5. H{a, q; a, 7, (5, z) is the 
second-order Fuchsian equation with four regular singu- 
lar points. One application of Heun function is to the 
Calogero-Moser-Sutherland System [29]. 

The Heun confluent function Hc{a, P,"f,6,ri, z) is ob- 
tained from Heun function through a confluence process. 
So Hc{a, P, 7, 6, 77, z) has two regular singular points and 
one irregular one. It is the singular points most impor- 
tant to our application. Hc(a,l3,^,5,rj,z) is defined as 
the solution of the following equations 

£■ d 
z{z - l)—r + [az^ + (-a -|-/3 + 7 + 2)z-/3 - 1] — -f 
rfz^ dz 

{a{p + 7 + 2) + 2d)z ~ a{(3 + 1) + /?(7 + 1) + 7 + 2?? 

2 

Hc{a,(3,'j,6,ri,z) = 0, 
Hc{a,f3,j,S,T],0) = 1, 

dH{a, 13, 7, (5, 77, z) Pi^a + + 1) - a + j + 2r] 



dz 



2(/?+l) 



(A2) 



The Heun doubly confluent function Hd{a, f3,"f, S, z) is 
obtained from Heun function through two confluence pro- 
cess and has two irregular singular points. Hd{a, (3, 7, 5, z) 
is defined as the solution of the following equations 



(= 



(2z^ 



+ 2z + a) — 



dz 



+ {f3z^ + {2a + 7)z + S) Ha{a, (3, 7, S, z) - 0, 



Hd{a,/3,j,S, z) 



dz 



z=0 = 0. 



(A3) 
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